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I.  INTRODUCTION 


Two  previous  DesmaCics  technical  reports  [3,  5]  discussed  impact 
acceleration  injury  prediction  model  development  and  estimation  accuracy 
respectively.  Another  technical  report  [4]  described  the  construction  of 
Impact  acceleration  injury  prediction  models  from  a set  of  twenty-eight 
-Gx  accelerator  rune  Involving  Rhesus  monkeys  with  securely  restrained 
torso  and  unrestrained  head.  The  models  described  in  that  report  employed 
head  dynamic  response  variables  and  sled  profile  variables  as  injury  predic- 
tors. 

This  report  considers  the  application  of  such  models  to  the  same 
twenty-eight  -C  acceleration  run  data  set  described  in  [4],  but  with  torque 
and  force  variables  as  possible  injury  predictors.  The  purpose  of  this 
report  is  to  determine  whether  the  forces  or  torques  experienced  during 
these  runs  can  be  used  to  predict  injury  probability. 

The  models  under  consideration  are  based  on  the  assumption  of  an  under- 
lying functional  relationship  of  the  form 

k -1 

P(x)  - (1  + exp[-(B0  + 1$^)]}  1 
where: 

x ■ (x^,...,x^)  denotes  the  set  of  independent  variables  considered, 
(Bq,...,^)  denotes  a set  of  parameter  values, 
and  P(x)  denotes  the  true  probability  of  injury  corresponding  to  x.  The 
torque  and  force  variables  used  in  this  investigation  were: 


(1)  peak  torque  around  the  anatomical  Y axis 


(2)  peak  fores  along  Che  anatomical  X axis, 
and  (3)  paak  force  along  the  anatomical  Z axis. 


Thasa  thraa  variables  vara  Investigated  In  a preliminary  study  conducted 
by  Lustlck  and  Shlmp  [1]. 

The  statistical  analysis  used  In  constructing  the  Impact  acceleration 
Injury  prediction  models  Indicates  that,  of  the  three  variables  considered, 
the  paak  force  along  the  anatomical  Z axis  correlates  most  highly  with 
Injury.  In  five  of  the  six  cases  where  Injury  occurred,  this  force  was 
76.3  kilograms  or  larger. 


II.  MODEL  CONSTRUCTION 


The  data  base  used  In  model  construction  consisted  of  twenty-eight 
-Gx  accelerator  runs  on  Rhesus  monkeys.  The  experimental  runs  are  the  same 
ones  as  In  technical  report  [4],  but  here  different  Impact  acceleration 
Injury  predictors  are  used.  Because  some  monkeys  were  run  more  than  once, 
dependence  exists  In  the  data.  It  will  be  assumed,  however,  that  the  effects 
of  dependence  are  small.  In  any  event,  the  dependence  In  the  data  will  be 
responsible  for  a slightly  conservative  model,  l.e.,  one  that  would  predict 
probabilities  of  Injury  that  are  biased  upward. 

Since  It  Is  difficult  to  define  Injury,  fatality  was  the  criterion  used 
In  development  of  the  models  discussed  In  this  report.  So  In  actualltv  these 
models  are  fatality  prediction  models.  The  complete  data  set  is  given  In 
Figure  1.  In  this  figure  a 1 represents  a fatal  run  and  a 0 represents  a 
nonfatal  run. 

It  should  be  noted  that  most  fatalities  Involved  a transection  in  the 
region  between  the  lower  medulla  and  upper  cervical  spinal  cord.  (See  (7] 
for  a further  discussion  of  the  neuropathologlcal  findings.)  This  la  where 
the  torque  and  force  variables  were  measured.  Thus,  the  models  in  this  report 
are  used  to  predict  fatality  from  torques  and/or  forces  at  the  site  of  the 
inj ury . 

A.  ESTIMATION  OF  MODEL  PARAMETERS 

A computer  program  for  maximum  likelihood  estimation  was  used  to 


calculate  S , 3 ,...,£,  l.e.,  the  estimates  of  the  parameters 

0 1k  0 i ^ 

In  model  construction,  it  la  possible  that  some  or  all  of  the  candidate 
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Figure  1:  The  Data  Set 
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variables  may  prove  unimportant  and  should  therefore  not  be  Included  In  a 
final  model.  The  contribution  of  variables  may  be  discerned  by  likelihood- 
ratio  tests.  These  tests  may  be  used  in  conjunction  with  "nested"  models. 

In  this  technical  report,  one  model  will  be  said  to  be  nested  within 
another  if  the  second  model  contains  all  variables  of  the  first  model  plus 
one  or  more  additional  variables.  Thus,  a model  which  contained  variables 
Xl’  *2'  *n<*  *3  wou^d  k®  nested  within  a model  which  contained  only  x^  and  x^. 

To  test  a hypothesis  that  a model  containing  variables  (x^ x^  + is  a 

statistically  significant  improvement  over  a model  containing  variables 

(x^ x^),  a log-likelihood  statistic  may  be  used.  The  procedure  Is  to 

calculate: 

• -2  log  likelihood  for  the  model  containing  (x^ x^) 

and  L ■ -2  log  likelihood  for  the  model  containing  (x, , ). 

2 1 K *r  tn 

Under  the  null  hypothesis  that  the  m additional  variables  (x^  + , ...,x^  + 

do  not  result  In  an  improved  model,  the  statistic  has  an  approximate 

Chi-square  distribution  with  m degrees  of  freedom.  Thus,  the  hypothesis 
may  be  tested  by  comparing  the  value  of  with  the  upper  percentage 

points  of  a Chi-square  distribution  with  m degrees  of  freedom. 

S.  TORQUE  AND  FORCE  VARIABLES 

As  mentioned  in  the  Introduction,  three  fatality  predictors  were  con- 
sidered for  model  development.  These  variables  will  be  denoted  by  THB,  FHA, 
and  FHC  (as  defined  by  Lustick  and  Shimp  [1]),  where: 

THB  is  the  peak  torque  around  the  anatomical  Y axis. 
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FHA  ia  the  peak  force  along  the  anatomical  X axis, 
and  FHC  la  the  peak  force  along  the  anatomical  1 axis. 

Model  construction  was  somewhat  complicated  by  two  aspects  of  the 
data  set.  First,  the  torque  variable,  THB,  for  observation  1365,  subject 
3921,  is  not  available.  Second,  Lustick  and  Shlmp  [1]  state  that  subject 
3935  may  have  been  killed  by  a combination  of  high  torque  along  the  Y axis 
and  high  force  along  the  X axis;  the  subject  being  compromised  on  run  1362 
and  finally  killed  on  run  1363.  This  suggests  that  observations  1362  and 
1363  may  not  be  homogeneous  with  the  rest  of  the  data  set. 

Therefore,  all  models  that  Include  the  THB  variable  were  estimated 
without  observation  1365.  Further,  all  models  were  estimated  with  and 
without  observations  1362  and  1363  to  observe  the  differences  in  the  result- 
ing models  and  statistics.  For  reference  purposes,  a data  set  that  contains 
observations  1362  and  1363  will  be  referred  to  as  an  A data  set;  otherwise 
the  data  set  will  be  referred  to  as  a B data  set.  For  a specified  number 
of  variables,  that  model  which  yielded  the  smallest  -2  log  likelihood  value 
was  selected  as  the  best.  As  can  be  seen  from  Figure  2,  the  best  one-variable, 
two-variable,  and  three-variable  models,  estimated  from  data  set  A,  are 
based  on  respectively, 

(1)  FHC 

(2)  THB,  FHC 

and  (3)  THB,  FHA,  FHC. 

Figure  3 implies  essentially  the  same  results  for  data  set  B.  Since  these 
models  are  nested,  the  relative  contribution  of  variables  FHC,  THB,  and  FHA 
may  be  tested  in  that  order. 
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Figure  2:  Head  Dynamic  Response  Variable  Sets  and  Associated 
-2  Log  Likelihood  Values  for  Data  Set  A 
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Figures  4 and  5 are  summaries  of  the  relevant  test  procedure.  In 
the  first  stage,  for  both  data  sets  A and  B,  FHC  was  tested  to  determine 
whether  it  significantly  Improved  the  model  which  assumed  a constant  prob- 
ability over  all  values  of  the  torque  and  force  variables.  Data  set  A 
yielded  an  observed  chi-square  statistic  of  18.46,  which  is  statistically 
significant  at  the  0.0000+  level.  Data  set  B yielded  an  observed  chi-square 
statistic  of  25.46  which  is  significant  at  the  0.0000+  level.  Thus,  both 
data  sets  indicate  that  FHC  significantly  contributes  to  fatality  prediction. 

The  second  stage  of  model  testing  considered  the  addition  of  another 
variable  to  the  model  which  Included  only  FHC.  Since  the  model  based  on  THB 
and  FHC  is  the  best  two  variable  model,  the  effect  of  adding  THB  to  the  model 
was  examined.  Including  the  THB  variable  resulted  in  observed  chi-square 
statistics  of  0.54  for  data  set  A and  0.00  for  data  set  B,  both  of  which  are 
not  statistically  significant.  Finally,  the  model  based  on  all  three  variables 
THB,  FHA,  and  FHC,  when  tested  against  the  model  based  on  FHC  only,  resulted 
in  statistically  insignificant  chi-square  statistics  of  0.56  for  data  set  A 
and  0.00  for  data  set  B. 

Thus,  based  on  both  data  sets  A and  B,  a prediction  model  which  includes 
only  the  FHC  force  variable  is  "best"  as  previously  defined.  The  resulting 
"best"  model,  estimated  from  data  set  A,  is  given  by 

- -1 

Model  A:  P(FHC)  - {l  + exp[-(-7.3705  - 0. 10721 (FHC) ) 1 } 

The  resulting  "best"  model,  estimated  from  data  set  B,  is  given  by: 

Model  B:  P(FHC)  - (l  + exp[-(-281.42  - 3.7902(FHC))]}_1 
Figures  6 and  7 show  a comparison  of  observed  probability  (i.e.,  0 or  1) 


Test  1:  FHC  against  Constant  Only 

^ - L2  - 29.10  - 10.64  - 18.46  (1  d.f.,  p - .0000+) 

Significant 

Test  2:  (THB,  FHC)  against  FHC 

L1  - L2  - 10.64  - 10.10  - 0.54  (1  d.f.,  p > .50) 

Nonsignificant 

Test  3:  (THB,  FHA,  FHC)  against  FHC 

- 1*2  ■ 10.64  - 10.08  = 0.56  (2  d.f.,  p > .50) 

Nonsignificant 

THB  denotes  peak  torque  around  the  anatomical  Y axis 
FHA  denotes  peak  force  along  the  anatomical  X axis 
FHC  denotes  peak  force  along  the  anatomical  Z axis 


Figure  4:  Testing  the  Significance  of  the  Head  Dynamic 
Response  Variables  Using  Data  Set  A 
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Teat  1:  FHC  against  Constant  Only 

4 * L2  - 25.46  - 0.00  - 25.46  (1  d.f.,  p - .0000'*') 

Significant 

Teat  2:  (TUB,  FHC)  against  FHC 

L1  “ L2  “ 0,00  " 0,00  * °’00  C 1 d-f*.  P > -50) 

Nonalgnlf leant 

Teat  3:  (THB,  FHA,  FHC)  agalnat  FHC 

Lj  “ L2  - 0.00  - 0.00  - 0.00  (2  d.f.,  p > .50) 

Nonsignificant 


THB  denotes  peak  torque  around  the  anatomical  Y axis 
FHA  denotes  peak  force  along  the  anatomical  X axis 
FHC  denotes  peak  force  along  the  anatomical  Z axis 


Figure  5:  Testing  the  Significance  of  the  Head  Dynamic 
Response  Variables  Using  Data  Set  B 
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Figure  7:  A Comparison  of  Observed  and  Predicted 

Probabilities  for  Model  B Based  on  Data  Set  B. 
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and  predicted  probability  for  data  seta  A and  B respectively.  In  both  of 
these  figures  the  observations  are  arranged  In  order  of  increasing  predicted 
probability. 


C.  PREDICTION  OF  CRITICAL  ENVELOPES 

Utilizing  models  A and  B,  a critical  envelope  for  each  can  be  predicted. 
This  envelope  defines  those  combinations  of  Independent  variables  for  which 
the  predicted  probability  of  injury  (or  fatality)  is  greater  than  some  speci- 
fied amount.  This  type  of  critical  envelope  corresponds  to  the  type  defined 
in  a previous  Desmatics  technical  report  [4]. 

Suppose  it  were  desired  to  restrict  the  variable  values  to  a region  in 
which  the  probability  of  fatality  were  less  than  some  small  probability  p 

A W 

(such  as  .01  or  .05).  In  other  words,  the  predicted  probability  £(x)  is 
desired  to  be  less  than  p^.  From  this,  the  three  following  equivalent 
inequalities  follow: 

(1)  P<£>  1 P0  , 

a k a ^ i 

(2)  (1  + exp[-(6Q  + Z 3^)]}  £ P0  , 

and  (3)  + E 5^  < ln[pQ/(l  - pQ)J  . 


For  model  A,  the  critical  envelope  at  p^  » .05  is 
-7.3705  - 0. 10721 (FHC)  < -2.9444 
or  equivalently 

FHC  > -41.29. 

For  model  B,  the  critical  envelope  at  p - .05  is 
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-281.42  - 3. 7902(FHC)  < -2.9444 
or  equivalently 

FHC  > -73.47  . 

Thus,  the  probability  of  fatality  is  predicted  to  be  less  than  52  if 
the  magnitude  of  the  force  along  the  negative  anatomical  Z axis  is  less 
than  41,29  and  73.47  kilograms  for  models  A and  B respectively.  Here  it 
is  easy  to  observe  that  the  exclusion  of  observations  No.  1362  and  No.  1363 
have  a large  effect  on  critical  envelope  prediction. 
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III.  SUMMARY  AND  DISCUSSION 

The  formal  statistical  analyses  of  the  data  on  the  twenty-eight 
Rhesus  monkeys  in  conjunction  with  previous  results  [4]  indicate  that, 
for  models  using  Rhesus  dynamic  response,  the  one  employing  only  the  FHC 
variable  (the  force  along  the  anatomical  Z axis)  is  the  best  model  for 
predicting  fatality  (injury).  However,  it  should  be  recognized  that  the 
total  amount  of  data  used  was  small  and  that  there  were  two  observations 
about  which  some  questions  exist. 

For  both  data  sets  A and  B,  the  "nested"  model  testing  indicates  that 
the  peak  torque  around  the  anatomical  Y axis,  THB,  and  the  force  along  the 
anatomical  X axis,  FHA,  do  not  contribute  significantly  to  fatality  predic- 
tion. The  testing  also  shows  that  the  FHC  variable  contains  the  most 
information  of  the  three  variables  under  consideration  whether  or  not  obser- 
vations No.  1362  and  No.  1363  are  included. 

For  observation  No.  1363  all  estimated  models  predict  low  probabilities 
of  fatality,  even  though  this  observation  corresponds  to  a fatality.  Although 
there  is  a possibility  that  this  observation  is  an  outlier  in  the  data  set, 
it  may  be  that  THB  and/or  FHA  are  also  important  predictor  variables,  but 
that  there  la  insufficient  Information  on  them,  l.e.,  not  enough  observations 
where  the  THB  and  FHA  values  are  high  in  the  negative  direction.  If  observa- 
tion No.  1363  is  Indeed  an  outlier,  then  the  "best"  model  estimated  from  data 
set  B,  model  B,  should  be  the  accepted  model  based  on  the  experimentation 
completed  thus  far.  However,  if  observation  No.  1363  is  a relatively  rare 
(but  Important)  observation,  then  model  A is  to  be  the  preferred,  parsimonious 
prediction  model. 
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Ac  present.  It  would  be  premature  and  unwise  to  conclude  that  observa- 
tions No.  1362  and  No.  1363  are  spurious  and  therefore  remove  them  from  the 
data  set.  It  should  be  noted  that  the  sled  profile  variables  provide  almost 
perfect  predictors  of  fatality  for  all  the  observations,  including  No.  1362 
and  No.  1363.  (See  (4].)  This  suggests  that  these  observations  should  not 
be  regarded  as  outliers.  Thus,  if  a model  were  to  be  used  at  present,  model 
A would  probably  be  the  wisest  choice.  It  should  be  beneficial,  however, 
to  attempt  to  extract  optimal  predictors  from  the  dynamic  response  by  the 
procedures  outlined  in  a forthcoming  report  [6], 
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